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I, INTRODUCTION 


Pneumatic tires are functional parts of many dynamic systems. In 
order to effectively design and engineer such systems, it is often neces- 
sary to know the mechanical properties of the component parts, including the 
tires. One of the major roles of this research group has been to study and 
analyze some of the important mechanical properties of pneumatic tires and 
to present rational methods for predicting them. 

Several paths have been followed by this research group in developing 
approximate methods for predicting mechanical properties of tires. One 
technique involved modeling of the penumatic tire as a cylindrical shell 

5-7 

supported by an elastic foundation. This model gives relations for pre- 

dicting several properties involving deformation in the plane of the wheel, 
such as contact patch length vs. vertical deflection, vertical load vs. 
vertical deflection, plane vibration characteristics, transmissability 
characteristics, and dynamic response to a point load. 

4 

A more recent attempt involved analyzing the pneumatic tire as a 
string on an elastic foundation. This model is primarily used to predict 
lateral stiffness characteristics, vertical stiffness characteristics and 
twisting moments. However, it has been found to be useful only when the 
inflation pressure is high, such as in the case of aircraft tires. 

This report presents a method for predicting the fore-and-aft stiffness 
characteristics of pneumatic tires. Fore-and-aft properties are important 


in the overall analysis of a tire since they represent the contributions of 



the carcass and tread to braking and tractive elasticity* A different model 


is required here since neither the cylindrical shell nor the string on the 
elastic foundation provide means for transmitting such loads* It is hoped 
that this model will prove satisfactory for predicting fore-and-aft charac- 
teristics of various tire designs* 
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II. SUMMARY 


An elastic bar supported by a foundation exhibiting elasticity in shear 
serves as a model for determining the fore-and-aft stiffness properties of a 
pneumatic tire. The differential equation representing the deformation is 
derived and solved, and the resulting solution gives a means for calculating 
a fore-and-aft spring rate for the model, 

A series of five tires of various sizes and structures was used for test- 
ing the validity of the proposed model, A set of static tests was performed 
to establish an experimental value for the fore-and-aft spring-rate for the 
various tires. Additional structural data, required by the analytical solu- 
tion of the model, were also obtained from the tires. A comparison of the 
calculated and experimental results was reasonably satisfactory, indicating 
that the proposed model can be used to roughly approximate fore-and-aft stiff- 
ness characteristics. 

A complete tabular summary of the geometry and composition of the five 
tires is included for easy reference. All experimental and analytical results 
are summarized and compared in graphical form. 
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III. ANALYSIS 


To represent fore-and-aft stiffness characteristics , the pneumatic tire 
is idealized as an elastic bar supported by an elastic shear foundation (see 
Figure l) . The elastic bar portion of the model represents the tread region 
of the tire which is loaded by the fore-and-aft load F. In addition to the 
restraint offered by the stiffness of the tread region itself, resistance to 
deformation by the load F is provided by the tires 1 ability to withstand 
shearing forces in the sidewall regions. This portion of the tire is repre- 
sented in the model by the elastic shear foundation. 



Figure 1. Idealized model of pneumatic tire for 
analyzing fore-and-aft spring rates. 

If it is assumed that the restraining force of the elastic shear founda- 
tion is directly proportional to the displacement, an elemental segment of 
the elastic bar can be set in equilibrium as shown in Figure 2. Note that ad- 
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vantage is taken of the symmetry present in the model 
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Figure 2. Loaded element of a portion of the model. 


In Figure 2 ; S is the force acting on the bar, u is the displacement, 
and K g is the spring rate per unit length of the shear foundation. It is 
assumed that K g is provided only by the shear resistance of the sidewall. 
From Figure 2, it is seen that one may approximate 

K 2GH 

K s - A 


where G is the effective shear modulus of the sidewall, H is the sidewall 
thickness, and A is the length along the sidewall from the rim to the point 
of intersection of the tread and carcass. 

From equilibrium of the element, 


dS 

3x 


K s u 


0 


S = TA S = EeA s = A S E 


5 


5u 

5x 




where T is the stress, e the resulting strain, A s the cross-sectional area of 


the bar at any location, and E the effective extension modulus of the tread 
region in the circumferential direction. Thus, 

~ - q. 2 u = o (i) 

ax 2 


where 


2 

<1 


K. 


s 



The general solution of this equation is 


u = C-^ cosh qx + sinh qx (2) 

The boundary conditions for this problem are determined by assuming that 
each half of the tire ( f ore- and- af t of the contact patch) is equally loaded, 
so that 

F 

at x = 0 , S = 2 

at x = jta, S = 0 (3) 


Substituting ( 3 ) into (2) gives 


-F 


2A g E tanh Ttaq 9 


°2 = 


2A s Eq 


Thus, 




(U) 


(5) 
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I 


The fore-and-aft stiffness is determined by finding the ratio of the applied 
load to the displacement at the point of application of the load, Thus, 


Kf 


F 


u(o) 



- 2A g Eq (tanh qira) 


(6) 


Equation (6) now represents a relationship for the fore-and-aft spring-rate 
of a pneumatic tire idealized as an elastic bar supported by a shear founda- 
tion,, As can be seen From Eqs„ (l) and (6), the application of Eq„ (6) to a 
real tire requires a knowledge of the effective stiffness (A s E) of the tread 
region in the circumferential direction, the effective shear modulus of the 
side-wall region G,the effective sidewall thickness H, and the length along 
the mean meridional section from the rim to the intersection of the tread and 
carcass, A 0 

The extension modulus in the circumferential direction and the shear 
modulus of the carcass usually vary from one location to another in the merid- 
ional direction because of the orthotropic nature of the tire carcass, so some 
criteria must be established to compute A g E and G for a given tire section 0 
Both the extension and shear modulus used for the results listed in this 
report were obtained by averaging the actual values of these properties 
throughout the cross-section,, This technique has been successful because 
the variation in these properties has not been too nonlinear* However, it 
has a great disadvantage in a simplified analysis such as this because it 
requires lengthy calculations which cannot be done efficiently without the 
aid of a digital computer,, For this reason an effort has been made to ob- 
tain some simplified, but reasonably accurate, approximations for the A g E 
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and G necessary for calculating the fore-and-aft spring constant. The re- 
sults of this effort are included in Appendix I and give a satisfactory ap- 
proximate technique for calculating these properties. 



-t»- X 


Figure 3. Meridional cross-section of the tire 
showing symbols used. 


In order to establish some validity for the assumption that the shear 
foundation modulus K„ can be estimated by considering shear effects only, a 

5 

simple experiment was performed by gluing a metal strip along the line of con 
tact of two rubber cylinders placed side by side (see Figure 4 ). A load was 
attached to the bar and the resulting deflection was measured by the dial in- 
dicator. The slope of the experimental load-deflection curve, related to K s , 
was then compared with the value of K s obtained from the relation given above 

v _ 2GH 
" A 


A summary of this experiment is presented below: 
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IV. COMPARISON OF THEORY WITH EXPERIMENT 


In order to investigate the validity of Eg.. (6), a series of static 
fore-and-aft stiffness tests were run on representative tires. Before 
reporting these tests and their results, the five tires used are described 
in detail. The idealized centerline profiles of the tires are shown in 
Figure 6. Tire No. 1 is a domestic 4-ply, 8.00 x 14 bias-ply tire with 
standard nylon cord. Tire No. 2 is a 2-ply, 7-50 x 14 bias-ply tire with 
standard nylon cord. Tire No. 3 is an imported 4-ply, 5*90 x 15 bias-ply 
tire with nylon cord. Tire No. 4 is an imported 7*50 x 14 radial-ply tire 
with overheads reinforced with wire cord. Tire No. 5 is a European made 
155 mm x 15 in. radial-ply tire with overheads reinforced with nylon cord. 

Table I is a summary of the pertinent elastic and geometric parameters 
required from the five tires. Using the results in this table, Figure 6, 

Eq. (6) and the proper elastic properties, it is possible to calculate the 
fore-and-aft stiffness of the five tires. Carrying out these computations 
gives the calculated values presented in Table II. 

To check the accuracy of the calculated values, the five tires were 
tested in the apparatus illustrated in Figures 7 and 8. In this testing 
procedure the tires were loaded vertically to a fixed deflection. Then a 
varying fore-and-aft load was applied and the resulting deflection recorded. 
The slope of these load-deflection curves represents the experimental fore- 
and-aft spring-rates. These tests were run for different vertical deflections 
and inflation pressures. The results of these tests are summarized in Fig- 
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Figure 6. Idealized mid-line profiles of five pneumatic tires. 










TABLE I 


SUMMARY OF GEOMETRIC, ELASTIC, AND STRUCTURAL PROPERTIES OF FIVE AUTOMOTIVE TIRES 




Item (Ref. 5, Table I) 

Tire 1 

Tire 2 

Tire 3 

Tire 4 

Tire 5 

Bias-Ply 

8.00x14 

Bias-Ply 

7-50x14 

Bias-Ply 

5*90x15 

Radial-Ply 

7.50x14 

Radial-Ply 
155mm xl5 

AO 

- 

outside radius of tire 

12.525 

13.9^ 

12.875 

13.44 

12.25 

L 

- 

half circumference 

59-55 

^3.79 

40.45 

42.22 

38.48 

ET 

- 

extension modulus, tread rubber 

670. 

560. 

481. 

690. 

490. 

A 

- 

length, mean meridional section 

4.6408 

5.84o4 

4.9217 

5.7238 

4.5529 

H 

- 

effective thickness for K s 

0.164 

0.110 

0.160 

0.250 

0.280 

G 

- 

effective shear modulus for K s 

28440. 

47164. 

43876. 

269. (G X y) 

144. (G xy ) 

a b e 

- 

effective spring rate, circumferential 

779. 

747. 

760. 

11906. 7 

23027. 


- 

spring rate, shear foundation 

2010. 

1747. 

2852. 

18.98 

23.50 

BETAC 

- 

cord half angle, crown 

0.6458 

0.6283 

0.6109 

0.3142 

0.2356 

ROC 

- 

radial location, crown 

12.120 

13.38 

12.355 

12.92 

11.70 

ROB 

- 

radial location, rim 

7.078 

7.20 

7.515 

7.06 

7.59 

R 

- 

idealized radius, sidewall 

5.22 

5-03 

2.84 

3-13 

2.17 

Yc 

- 

y-coordinate, center for R 

1.955 

2.490 

2.092 

2.454 

1.939 

x c 

- 

x-coordinate, center for R 

-0.057 

0.465 

-0.146 

0.405 

0.704 

R1 

- 

idealized radius, crown region 

5.40 

2.62 

2.76 

3-54 

4.88 

ERUB 

- 

extension modulus, carcass rubber 

438. 

310. 

370. 

625. 

300. 

GRUB 

- 

shear modulus, carcass rubber 

146. 

103- 

123. 

208. 

100. 

MURUB 

- 

Poisson ratio, carcass rubber 

0.500 

0.500 

0.500 

0.500 

0.500 

GCORD 

- 

shear modulus, cord 

705. 

705. 

705. 

705. 

705- 

AESUBC 

- 

spring rate, cord 

200. 

623. 

317. 

350. 

307. 

DIAMC 

- 

effective diameter, cord 

0.025 

o.oUo 

0.026 

0.025 

0.023 

MUC 

- 

Poisson ratio, cord 

0.700 

0.700 

0.700 

0.700 

0.700 

TPLY 

- 

effective ply thickness 

0.041 

0.055 

0.040 

0.040 

0.040 

NCORD 

- 

cord count, crown 

26. 

19. 

24. 

18. 

20. 

ALPSTR 

- 

normal angle, intersection 

0 . 7746 

0 . 6085 

0 . 6665 

0 . 6427 

0.5775 

ALPHR 

- 

normal angle, rim 

2.2148 

2.5551 

2.3973 

2.4714 

2.6792 

BW 

- 

tread width 

4.36 

4.40 

3.20 

4.60 

3.68 


TABLE II 

SUMMARY OF EXPERIMENTAL AND CALCULATED VALUES OF K f 


Tire 

Kf - 

lb/in. 

Experimental 

Calculated 

1 

2530 

2503 

2 

2300 

2286 

3 

2780 

2944 

4 

1340 

1469 

5 

1265 

950 
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ures 9 through 13 , 


In general it can be seen from these curves that the fore-and-aft spring- 
rate increases only slightly with increasing vertical load and with increasing 
internal pressure. Since Eq, (6) does not account for the slight increase 
due to these factors y a comparison between the experimental and calculated re- 
sults must be made in a somewhat arbitrary fashion. However, since the experi- 
mental values are nearly the same for all conditions examined, any results used 
as a comparison with the calculated values will serve as a meaningful check. 

The comparisons shown in Table II are based on experimental values obtained 
from vertical tire deflections of one inch and by use of manufacturers .rated 
inflation pressure. The experimental values were determined by measuring the 
slopes in the linear portions of the load-deflection curves. These compari- 
sons indicate that the simple model formulated above gives a method for ap- 
proximating the fore-and-aft spring-rate of pneumatic tires using only the 
geometric, elastic, and structural properties required by most tire designers. 
At the end of the Appendix II an example problem is worked out, illustrating 
how an approximate fore-and-aft spring rate can be calculated if the correct 
input data is available. 
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TRACTIVE FORCE- LBS 


o 1/2 in. vertical deflection 

x 3/2 in. vertical deflection 

a I in. vertical deflection 

□ 11/4 in. vertical deflection 



Figure 10. Fore-and-aft load-deflection curves. Tire No. 2. 


o I /2 in. vertical deflection 
x 3/4 in. vertical deflection 

a I in. vertical deflection 

□ 11/4 in. vertical deflection 



Figure 11. Fore-and-aft load-deflection curves. Tire No. 3. 


o 1/2 in. vertical deflection 

x 3/4 in. vertical deflection 

a I in. vertical deflection 

□ I 1/4 in. vertical deflection 



HORIZONTAL DEFLECTION - IN. 


Figure 12. Fore-and-aft load-deflection curves. Tire No. 4. 



il deflection 
>1 deflection 
il deflection 
I deflection 



Tire No. 5 



V. APPENDIX I 


This Appendix presents the results of an effort to simplify the compu- 
tation's involved in obtaining the extension modulus of the carcass in the 
circumferential direction, Eg, and the shear modulus of the carcass, G. 

Reference 1 gives exact expressions for the moduli of laminated ortho- 
tropic two-dimensional sheets, and these are good representations for the 
elastic constants of a tire carcass. However the expressions derived in 
Ref. 1 are quite lengthy to evaluate. They can be simplified considerably 
by considering the structure to be made of inextensible cords, so that the 
modulus of elasticity E^, parallel to the cords in a single ply, becomes 
indefinitely large. By simplifying the expressions of Ref. 1 in this way, 
one gets 

UG VV sin 2 a cos 2 a + E„(cos 2 a - sin^a) 2 

E e - (7) 

sin a 

where a is the local cord half-angle at any meridional location, G is the 
shear modulus of an individual ply and is the extension modulus of an in- 
dividual ply in the direction normal to the cords. This is a very gcod ap- 
proximation for Eg as long as the cord half-angle is not less than 30 °. This 
is shown in Figure 14 where Eq. (7) is compared with the exact formulation of 
Ref. 1. Thus, Eq. (7) is a good approximation for an ordinary bias-construc- 
tion tire, where the angle is almost always greater than 30 ° • 

An approximation for the shear modulus can be obtained from a strength 
of materials analysis. This is reproduced in Appendix II in detail, and gives 
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Circumferential Extension Modulus E 






G 


(8) 


2 2 

= E^sin & cos a 


where E^. is the extension modulus in the direction of the cords of an 
individual ply of the tire carcass material. This is a very good approxima- 
tion for the shear modulus except at cord half-angles near 0 and 90° as 
shown in Figure 15 where the exact shear modulus expression from Ref. 1 
is compared with Eq_. (8). Again, this is valid approximation for almost 
all tire constructions since cord angles usually lie between 30° and 60°. 

Reference 2 presents a concise method for calculating E^, Ey, 

and G as a function of the geometric and elastic properties of a single 
xy 

ply of the carcass material. The expressions for these moduli are listed 
below : 


E x = ( AESUBC) ( NCORD) /TELY 




= ERUB 


1 A \ 

1 + 2 -9(-^ l y 5 ' 


G xy = 703 K d \ H \ s + GHUB(l-\ s X H ) 


where : 


K( DIAMC) 

= " W ” 


\ s = K( DIAMC) (NCORD) 


The constant k is an area coefficient equal to s/ it Jb. The other 
parameters are defined in Table I of this report. Unfortunately both the 
expression for E and G are still functions of the cord half-angle, a, 
which varies from position to position around the cross-section. However, 
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it has been found possible to select average Locations at which the values 


of the extension modulus and shear modulus represent useful values for the 
entire cross-section. 

The extensional stiffness of the "bar" is represented by A g E and 
primarily depends on the extensional stiffness of the carcass in the tread 
shoulder region of the tire D The cord half -angle in the tread and shoulder 
region of the tire ordinarily varies from about 35° to 45° and since the 
computation for is very simple at 45° , an effective E^ is selected as 
the value corresponding to the value of Eq at 45°. The total A S E is then 
calculated by multiplying the E Q at 45° by the total carcass thickness and 
the width of the tread. Thus,, 

.A S E - Eq(45°) • H • BW (LO) 

The shear modulus of the supporting foundation, G, is the shear 
modulus of the sidewall portion of the tire. Therefore G is approximated 
by its value at a location mid-way between the rim and the shoulder. in 
Figure 16 this position is noted by the dimension FS. The cord-half -ang Le , 
a s , is approximated at this position by the "cosine law": 

pq 

cosa s = r5c cosB ( 1L ) 

Where £3 is the cord half-angle at the crown. 

The relative accuracy of these a.pproximate moduli is examined by 
comparing their values with those used previously for Tires 1, 2, and 3 
in Table 1. These comparisons and the resulting effects on the fore-and- 
aft spring constant, are shown in Table IU- 
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TABLE III 


COMPARISON OF RESULTS USING APPROXIMATE MODULI 



Tire 1 

Tire 2 

Tire 3 

G'— From Table I 

28440 

47164 

43876 

G— Approx, from Eqs. (8) and (ll) 

29000 

48500 

45400 

A g E — From Table I 

779 

747 

760 

A S E — Approx, from Eq. (7) - a = 45° 

780 

580 

533 

Kp — From Table II 

2503 

2286 

2944 

— Using Approx. A E and G 

2530 

2070 

2520 

K-p — Experimental 

2530 

2300 

2780 




As can be seen, the comparisons are relatively good and since this 
model for calculating is simple, it seems justifiable to use these 
simpler expressions for the moduli in approximate calculations of the 
spring rate K^.. 

To summarize the method for calculating and to illustrate the 
relative ease with which it can be done if one uses the approximate moduli 
outlined above, an example problem is given below. 


PROBLEM STATEMENT 


Calculate an approximate fore-and-aft spring rate for the tire described 
in Figure 17 . This tire roughly corresponds to a standard 9*50 by 14 
four-ply tire. 

The elastic properties of the individual ply are calculated first 
by referring to Eqs. (9): 


X H = -T355T = 0-75B 


>- B =/£(• 028 ) (20) = 0.500 


E = ( 300 ) (20)/( . 053 ) = 1.82 x 1CT psi 


=sr = 


= 550 


1 + 2.9 


/ 0. 758 x 0. 500 \ 
V 1 - 0.500 ) 


= 1760 psi 


= 705 ( 0 . 785 ^) (0.758) (0.500) + 180 (1-0.758 x 0.500) 


= 325 psi 


The extension modulus E Q is then computed from Eq. (7)> using a = 45°: 


27 



AESUBC =300 lb/ in. 
ERUB = 550 psi 
BW = 5.35 in. 

ROC = 14.25 in. 
DIAMC =0.028 in. 


NCORD =20 ends/ in. 
GRUB = 180 psi 
p = 37.5° 

RI- 13.30 in. 
a = 14.83 in. 


TPLY =0.033 in. 
H = 0.132 in. 
ROB = 7.30 in. 

A = 9.70 in. 


Figure 17 . Description of hypothetical tire for example problem. 



4(325) 


1300 psi 




Thus A E can be found from Eq. (10) 


A S E = 1300 ( 0 . 132 ) (5-35) - 915 lb 


The next property to calculate is G, whose mean value is found at a 
location approximately half way between the rim and the shoulder. The 
radial location RS is determined by locating the point half way between 
ROB and PI (see Figure 16) , which for this example is: 


RS =• 


RI +■ ROB 


7-30 + 13-3 0 
2 


10.30 in. 


The cord half-angle at this location is found from Eq. (ll): 


cos a s - cos 37.5° = 0.573 


and sin a_ = 0.819. Thus from Eq. (8): 

o 


G = 1.82 x 10 5 (0,573) 2 (0.8l9) 2 = 4.00 x 10^ psi 


Now referring to Eq, (6), the fore-and-aft spring rate is 


K f = 2( A g E) (q) ( tanh qna) 


The value for q is found from Eq. (l) and the definition of K 


Ks „ 2GH 2(4.00 x 10 ) (0.132) = 1070 lb/in „ /jLn . 


q. = " 1 . 08 '/in. 
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Therefore, the spring rate is 


K f = 2(915) (1. 08) (1) 


I 98 O lb/in. 
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VI. APPENDIX II 


Consider a plane element made up of a series of parallel cords, with 
end count n, each carrying a tension load T Q as shown below. 



The number of cords per unit length 
on the upper and lower faces is n * 
cos a. The tension component per 
cord tangential to the upper or 
lower faces is T 0 sin a . Hence, 
the shear force per unit length 
is 


T 0 sin a • n cos a. 


If the thickness of the lamina is h, then the shear stress t, which by 
equilibrium acts on all edges equally, is 

t = t q tj sin a cos a (12) 

Now consider deformation of the element shown above in the shear direc- 
tion, based on the concept that the cord will elongate under tension T 0 . 
From the geometry. 


x sin a = 

where € is cord strain, i Q the original cord length. 

x = €£ 0 /sln a . 
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But shear strain is defined as 


6 i. 


^ i D cos a 


i Q sin a cos a 


sin a cos a 


(13) 


Also, cord tension and strain are related by the cord spring rate 


T 0 = (AE) C . € 

Hence, shear modulus G is 

T ^ sin a cos a 

G = — = / (AE) sin a cos a 

7 T 0 c 

= ( AE : ) c -■ sin 2 a cos 2 a (14) 

But the extension modulus parallel to the cords is (AE) C “ = E x , so 
that 

p o 

G = E x sin a cos a ( 15 ) 
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